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Abstract 

The famous F5 algorithm for computing Grobner basis was presented by Faugere in 2002. 
The original version of F5 is given in programming codes, so it is a bit difficult to understand. 
In this paper, the F5 algorithm is simplified as F5B in a Buchberger's style such that it is 
easy to understand and implement. In order to describe F5B, we introduce F5-reduction, 
which keeps the signature of labeled polynomials unchanged after reduction. The equivalence 
between F5 and F5B is also shown. At last, some versions of the F5 algorithm are illustrated. 
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1. Introduction 

Solving systems of polynomial equations is a basic problem in computer algebra, through 
which many practical problems can be solved easily. Among all the methods for solving 
polynomial systems, the Grobner basis method is one of the most effi cient approaches. After 



the concept of Grobner basis is proposed in 1965 (IBuchb 



erger 



been presented for compu t ing Grobner basis , including (ILazardl. Il983l ; iGebauer and Moller 



19861 : iGiovini et all Il99lb Mora et all 11992 : iFaugereL Il999l |2002|). Currently, F5 is one of 



1965), many algorithms have 



the most efficient algorithms. 

After the F5 algorithm is proposed, many researches have been don e. 
Bardet et al. study the complexity of this algorithm in (IBardet et all 12 004). 



For example, 
F auger e 



and Ars use the F5 algorithm to attack multivariable s ystems i n (IFaugere and Arsl . 120031 ). 
Stegers revisits the F5 algorithm in his master t hesis (jStegersl . 120051 ). Eder discusses the 
two crite ria of the F5 algorithm in ( lEderl . 120081 1 and proposes a variation of the F5 al- 
gorithm (lEder and Perrvl . l2009i ). Ars and Hashemi present two variation of c riteria in 



(lArs and Hashemi! . 120091 ) . Recently, Gao et al. give a new incremental algorithm in ( )Gao et al. 



20101 ) . The cu rrent authors discuss th e F5 algorithm over boolean ring and present a branch 
F5 algorithm (ISun and Wang! . l2009allbh. We also give a complete proof for the correctness 
of the F5 or F5-like algorithm in fjSun and Wang . 120101 ) . 
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However, since the original F5 algorithm is reported in programming codes, it seems a 
bit difficult to understand this famous algorithm very well. In this paper, we revisit the F5 
algorithm from our perspective and we simplify F5 to F5B in Buchberger's style, which is 
equivalent to the original F5 algorit hm. We also d iscuss some versions of the F5 algorithm, 
such as IncrementalF5 algorithm in (jFaugerd . 120021 ). Non-incremental F 5 algorithm reporte d 
by Faugere in INSCRYPT 2008 and Matrix F5 algorithm mentioned in (IBardet et all I2004J ). 

This paper is organized as follows. The main idea of the F5 algorithm is illustrated in 
Section 3 after some basic notions given in Section 2. In Section 4, F5B is proposed and the 
equivalence of F5 and F5B is shown as well. Some versions of the F5 algorithm are discussed 
in Section 5. At last, some conclusions are presented in Section 6. 



2. Basic Notations 

Let K be a field and K.[X] = K[x\, ■ ■ ■ , x n ] a polynomial ring with coefficients in K. Let 
N be the set of non- negative integers and PP(X) the set of power products of {xi, ■ ■ • , x n }, 
i.e. PP{X) := {x a | x a = x^ 1 ■■■x^,a i 6 N, % = 1, • • • , n}. 

Let -< be an admissible order defined over PP(X). Given t = x a G PP(X), the degree 
of t is defined as deg(i) := \a\ = Y^i=i ai - For a polynomial 0^/6 K[xi, • • • , x n ], we have 
/ = Yl c a x a . The degree of / is defined as deg(/) := max{|a|, c a ^ 0} and the leading power 
product of / is lpp(/) := max^{x Q , c a ^ 0}. If lpp(/) = x a , then the leading coefficient and 
leading monomial of / are defined as lc(/) := c a and lm(/) := c a x a respectively. 



3. Revisit the F5 Algorithm 

In brief, the major contribution of the F5 algorithm is presenting two new criteria: 
Syzygy Criterion and Rewritten Criterion. Syzygy Criterion is also called F5 Criterion in 
some papers. Almost all useless computations (redundant S-polynomials) can be removed 
by these two criteria. Both of the criteria are built on the concept of signature and a special 
reduction procedure. 

In this section, we first illustrate the main idea of Syzygy Criterion and Rewritten Cri- 
terion by some simple examples, and then introduce the definition of signature. The special 
reduction is pre sented in next section. For more details about the F5 algorithm, please see 



flFaugerd . 120021 ) 



3.1. About Syzygy Criterion 

Of the two criteria, Syzygy Criterion is more important and creative. Let us see a simple 
example. 

Example 3.1. Consider a system in Q[x,y,z] with the Graded Reverse Lex Order with 

x y y y z, 

fi = x 2 + y, 
h = xy- z. 
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The S-polynomial spoly(/i, f 2 ) = yfi — xf 2 = y 2 + xz, which is not reducible by {/i, / 2 }. 
We denote / 3 = spoly(/i, /2). Now, consider the S-polynomial of / 3 and / 2 : 



the leading monomial of / 2 and / 3 are not co-prime, but spoly(/ 3 , / 2 ) reduces to by {fi, f 2 }. 
Why? What makes this S-polynomial spoly(/ 3 , / 2 ) reduce to 0? 

If we dig it deeper, we will see that the reason still comes from the syzygy (/ 2 , —fi) of 
the 2-tuple vector (/i, / 2 ). Since / 3 = y/i — x/ 2 and f 2 fi — /i/ 2 = 0, then we have 

spoly(/ 2 , /s) = xf 3 - yf 2 = xyf\ - (x 2 + y)/ 2 = lm(/ 2 )/i - /1/2 = (M/2) ~ /2)/i- 

Thus, it is natural that the S-polynomial spoly(/ 3 , / 2 ) reduces to by {/1}. 

Theoretically, in order to speed up the algorithm for computing Grobner basis, we have 
to avoid computing the S-polynomials of the above kind. Now, the question is: how to 
detect them? As the S-polynomial spoly(/ 3 , / 2 ) reduces to due to the syzygy (/ 2 , —fi) 
of the vector (f\, f 2 ), a natural idea is to connect the polynomials / 3 and f±. 

The relation / 3 = yfi — xf 2 means / 3 comes from yf\, so we can append a signature 
yei to / 3 to reflect this fact. Similarly, we also can append a signature xye\ to xf 3 , which 
shows xfz comes from xyf\ = lm(/ 2 )/i. Now from the signature xye\ = lm(/ 2 )ei, we are 
able to understand why the S-polynomial spoly(/ 3 , f 2 ) can be reduced to 0. 

In more general cases, for the polynomial system {/1, • • • , / m } C K[X], the syzygies of 
the m-tuple vector • • • , f m ) e (K[X]) m also result in many S-polynomials reducing to 0. 
So the Syzygy Criterion is a criterion that detects useless S-polynomials by using syzygies of 
(fi, • • • , fm)- In the above simple example, the signature of polynomial provides an useful 
information to detect this kind of unnecessary S-polynomials. Next, we give a mathematical 
explanation of signatures. 

3.2. Signatures and Labeled Polynomials 

Consider a polynomial system {/1, • • • , f m } C K[X] and denote (/1, • • • , /m) a polynomial 
m-tuple in (K[X]) m . We call the /j's initial polynomials of the ideal {fi, • • • , f m ), since they 
are initial generators of ideal • • • , / m ) C K[X]. 

Let ej be the canonical i-th unit vector in the free K[X]-module (K[X]) m , i.e. the i-th 
element of e^ is 1, while the others are 0. Consider the homomorphism map a over the free 



Then c(ej) = fi. More generally, if g = g\e\ + • • • + g m e mi where gi G K[X] for 1 < i < m, 

then cr(g) = yi/i H h g m f m . 

The admissible order -< on PP(X) extends to the free module (K[X]) m naturally in a 
POT (position over term) fashion: 



spoly(/ 3 , fz) = xf 3 - yf 2 = x 2 z + yz = zj u 



K[X]-module (K[X]) m : 



a : ( K [X]r — »- K[X], 

(^l> ' ' ' > fi 1 ™) 1 >■ ^l/l H 1" 9mfm- 
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Thus we have e m < e m _i < • • • < ej directly. 

With the admissible order < on (K[X]) m , we can define the leading power product, 
leading coefficient and leading monomial of a m-tuple vector g G (KLY]) m . For example, 
let g = (2x 2 + y 2 , 3xy) G (Q[x, y]) 2 or equivalently g = (2x 2 + y 2 )e\ + 3xye 2 and -< the Lex 
order on PP(x,y) (x y y). Then we have lpp(g) = x 2 e 1; lc(g) = 2 and lm(g) = 2x 2 B\. 

The following is the definition of la beled polynomia l and its signature. The idea of 



labeled polynomials draws lessons from fjGao et all 120101 ) 



Definition 3.2 (labeled polynomial). Let g G (fi, - ■ ■ , f m ) be a polynomial and g G (K[X]) m 
an m-tuple vector such that er(g) = g. Then we call Q = (g,g) a labeled polynomial. 
For a labeled polynomial Q , we define 

1. the signature sign(^) := lpp(g) 

2. the polynomial part poly(^) := g 

3. the leading power product lpp({?) := lpp(g) 

4. the leading monomial lm(Q) := lm(g) 

Suppose J 7 = (f , /), Q = (g, g) are labeled polynomials and cx a is a non-zero monomial, 
we define scaler multiplication and addition for labeled polynomials as following. 



• cx a ■ T = cx a T = (cx a f , cx a f). 

• F + G = (f + &,f + g). 



Definition 3.3 (critical pair, S-polynomial). For labeled polynomials J 7 , Q, we say [J 7 , Q] := 
(u, J 7 , v, Q) is a critical pair of J 7 and Q ifu,v are monomials such that ulm^J 7 ) = t>lm(C?) = 
lcm(lpp(J r ), lpp(C?)) ; and the corresponding S-polynomial of [J 7 , Q] is denoted by spoly(J r , Q) = 
uT — vQ . 

3.3. Syzygy Criterion 

We are now able to describe the Syzygy Criterion in mathematical words. We begin by 
the following definition. 

Definition 3.4 (Divisible). Let J 7 = (f , /) be a labeled polynomial with lpp(f) = x a ei, cx 1 a 
non-zero monomial and B a set of labeled polynomials. The labeled polynomial ex 1 J 7 is said 
to be divisible by B, if there exists a labeled polynomial Q — (g,g) G B with lpp(g) = x^ej 
such that 

1. Ipp(s0 | x 1+a , and 

2. e^ > Bj, i.e. i < j . 

Then the Syzygy Criterion is described as follows. 

Syzygy criterion is also called F5 criterion in some papers. In F5, if a critical pair meets 
the Syzygy Criterion, then the corresponding S-polynomial is redundant. 

Let us review the two simple examples from Subsection 13.11 in the view of Syzygy Crite- 
rion. 
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Criterion 1 — Syzygy Criterion 



Let [J 7 , Q\ := (u, J 7 , v, Q) be a critical pair and B a set of labeled polynomials. If either uT 
or vQ is divisible by B, then the critical pair [J 7 , Q\ meets the Syzygy Criterion. 



In Example 13.11 Let T\ = (ei,/i) and T 2 = (e 2 ,/2). Let J" 3 = spoly(J r i, J" 2 ) = y^i — 
2^2 = (2/ei - xe 2 , y 2 + xz). 

Next, consider the critical pair [^3,^2} = {x,3~3,y,J~2)- However, = (xyei — 

x 2 e 2 ,x(y 2 + xz)) is divisible by {J-2}, so the critical pair [^3,^2} is also redundant by 
Syzygy Criterion. And spoly(J-3, JF2) does reduce to in Example 13. 1[ which verifies the 
Syzygy Criterion. 

3.4- Rewritten Criterion 

We next describe the Rewritten Criterion. In fact, Rewritten Criterion is more like a 
programming technique, which reflects Buchberger's idea: try to reuse as much as possible 
the previous computations. Let us see the following example. 

Example 3.5. Consider a system in Q[x, y] with the Graded Reverse Lex Order with x y y: 

f fi = x 2 + xy, 
1 f2 = x 2 + y. 

The S-polynomial spoly(/i, / 2 ) = fi — f2 = xy — y is not reducible by {/1, / 2 }. So denote 
f'3 = xy — y. Next we need to consider the S-polynomial spoly(/i, fa) = yfi — xf% as well as 
the S-polynomial spoly(/ 2 , / 3 ) = yf 2 — xf 3 . However, with a further study, we will find that 
the S-polynomial spoly(/i, fa) is not necessary to reduce due to the following observation: 

spoly(/i, f 3 ) = yfi - xf 3 = y(fi - f 2 ) - (yf 2 - xf 3 ) = yspo\j(f u f 2 ) - spoly(/ 2 , / 3 ). 

That is, if the S-polynomials spoly(/i, / 2 ) and spoly(/ 2 , / 3 ) have been reduced (or scheduled 
to reduce), then reducing the S-polynomial spoly(/i, / 3 ) becomes a duplicated jobJl 

Thus, the Rewritten Criterion is a criterion that detects duplicated reductions of poly- 
nomials. The meaning of Rewritten Criterion is much clearer in a variation of the F5 
algorithm that uses matrix reduction, which will be discussed in Section |5j Now we give the 
mathematical definition of Rewritten Criterion. 

Definition 3.6 (Rewritable). Let T = (f , /) be a labeled polynomial, cx 1 a non-zero mono- 
mial in X and B a set of labeled polynomials. The labeled polynomial cx 1 ? = (cx 7 f , cx 1 f) 
is said to be rewritable by B, if there exists a labeled polynomial Q = {g,g) G B, such that: 

1. sign(^) I sign(ca; 7 J 7 ) , i.e. lpp(g) | lpp(cx 7 f), and 

2. labeled polynomial Q is generated later than T . 



5 



Criterion 2 — Rewritten Criterion 



Let [J 7 , Q] := (u, J 7 , v, Q) be a critical pair where u, v are monomials in X such that Mlm(J r ) = 
flm(^) = lcm(lpp(J r ), lpp(^)), and B a set of labeled polynomials. If either uT or vQ is 
rewritable by B, then the critical pair [J 7 , (?] meets the Rewritten Criterion. 



The Rewritten Criterion is given as follows. 

If a critical pair meets the Rewritten Criterion in the F5 algorithm, then it is not 
necessary to reduce its S-polynomial. Now let us explain the fact in Example 13.51 The 
system is fi = x 2 + xy,f2 = x 2 + y G Q[a:,y]. Labeled polynomials T\ = fi) and 
T 2 = (625/2) correspond to fi and f 2 respectively. The S-polynomial of T\ and T 2 is 
Tz = T\ - T 2 = (ei - e 2 , xy-y). 

Next, let us see the critical pair [Ti, J-3] = (y, T\, x, J-3). Clearly, labeled polynomial J-3 
is generated later than T\, so it is easy to know yT\ is rewritable by \T%\. Thus, critical 
pair [J^i, J-3] meets the Rewritten Criterion, and hence, it can be removed. 

As discussed in this section, both Syzygy Criterion and Rewritten Criterion build on the 
concept of signature. However, to ensure both criteria correct during the computations, a 
special reduction procedure is necessary, which is detailed in the next section. 



4. Simplify the F5 algorithm to a Buchberger's Style 

In order to simplify the F5 algorithm to a Buchberger's style, the most important step is 
to rewrite the reduction procedure. The original reduction in the F5 algorithm is described 
by codes and may return several reduction results each time. In this section, a simplified 
version of reduction (F5-reduction) is proposed in the first subsection. The F5 algorithm in 
Buchberger's style (F5B algorithm) is described in the second subsection. We will discuss 
the equivalence of the F5 and F5B algorithms in the last subsection. 

4-1. F5-Reduction 

The signatures of the labeled polynomials are used to detect useless critical pairs by the 
criteria, but this is not sufficient to ensure the correctness of the F5 algorithm. Only under 
F5-reduction, which is a special kind of reduction process, the critical pairs detected by the 
criteria are really useless. The same is true for other F5-like algorithms. 

Let us start with the definition of F5-reduction, which is a revised version of the TopRe- 



duction function in ( iFaugerd . 120021 ) . 



Definition 4.1 (F5- reduction). Let J 7 = (f , /) be a labeled polynomial and B a set of labeled 
polynomials. The labeled polynomial T is F5-reducible by B, if there exists Q = (g, g) G B 
such that: H 



2 Ri gorous proof needs th e knowledge of ^-representation. For more details, please see the authors' another 
paper (|Sun and WaneLl201Ch . 

3 Deleting the conditions 3 and 4 does not affect the correctness of algorithm, but leads to redundant 
computations/reductions. 
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1. lpp(fl-) I lpp(/), denote x 7 = lpp(/)/lpp(<7) and c = \c(f)/lc(g), 

2. sign(J r ) > sign(cx 7 (?) ; i.e. lpp(f) > lpp(cx 7 g), 

3. cx^Q is not divisible by B, and 

4. cx^Q is not rewritable by B. 

If T is F5-reducible by B, let J 7 ' = T — cx J Q. Then this procedure: J 7 J~' is called 

one step F5-reduction. If J 7 ' is still F5-reducible by B, then repeat this step until J 7 ' is not 
F5-reducible by B. Suppose J 7 * is the final result that is not F5-reducible by B. We say T 



F5-reduces to J 7 * by B, and denote it as T 



B 



J 7 *. 



The key of F5-reduction is the condition sign(J r ) > sign(cx 7 (?), i.e. lpp(f) >lpp(cx 7 g), 
which makes F5-reduction much different from other general reductions. The major function 
of this condition is to preserve the signature of J 7 during reductions. Thus a direct result is 
that, if labeled polynomial T F5-reduces to J 7 * by B (i.e. T J 7 *), then the signatures 
of J 7 and J 7 * are identical, i.e. 

sign (J 7 ) = sign (J 7 *). 

This property plays a cr ucial role in the proof for the correctness of the F5B algorithm. For 
more details, please see fjSun and Wang . I2010I ). 



4-2. The F5 algorithm in Buchberger's style 

With the definitions of Syzygy Criterion, Rewritten Criterion and F5-reduction, we can 
rewrite the F5 algorithm in Buchberger's style (F5B algorithm). 

Of course, some auxiliary data are also necessary to be recorded in the implementation. 
For example, we need to keep the generating order of the labeled polynomials. We can also 
save labeled polynomials as (lpp(g),5 f ) instead of (g,g) during the compuation, since only 
the information of lpp(g) is really used. 

The strategy of selecting critical pairs is not specified in the F5B algorithm, instead we 
simply use 

cpi — select a critical pair from CP. 

Different strategies of selecting critical pairs will lead to different versio ns of the F5 algo- 
rithm, including the original Incremental F5 algorithm in f Faugerel . 20021 ) . N on- Incremental 
F5 algorithm report ed by Faugere in re cent conference (INSCRYPT 2008), Matrix F5 algo- 
rithm mentioned in ( iBardet et al.l . 120041 ) . All these versions of the F5 algorithm are discussed 
in the next section. 

Moveover, the original F5 algorithm differs from the F5B algorithm only by a strategy 
of selecting critical pairs as well. Let us see more discussions in the following subsection. 

4-3. Equivalence of the F5 and F5B algorithms 



The major difference between F5 (jFaugerd . 120021 ) and F5B is the reduction process. In 
( IFaugerd . 120021 ). reduction of S-polynomials is done by the Reduction function. So next, we 
will focus on showing that the Reduction function is equivalent to a set of F5-reduction with 
an appropriate strategy of selecting critical pairs. Let us see the TopReduction function first, 
which is a subfunction of the Reduction function (function [1]). 
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Algorithm 1 — The F5 algorithm in Buchberger's style (F5B algorithm) 

Input: a polynomial m-tuple: (fi, ■ ■ ■ , fm) G -fT[X] m ; and an admissible order -<. 
Output: The Grobner basis of the ideal • • • , f m ) C 

begin 

Fi< — (e i( fi) for % = 1, • • • , m 

Bi — {Fi | i = 1, • • • ,m} 

CPi — {critical pair [J-^J-j] | 1 < i < j < m} 

while CP is not empty do 

cp< — select a critical pair from CP 

CP^CP \ {cp} 

if cp meets neither Syzygy Criterion nor Rewritten Criterion, 
then 

SP< — the S-polynomial of critical pair cp 
Vi — the F5-reduction result of SP by B, i.e. SP V 
if the polynomial of V is not 0, i.e. poly(P) ^ 0, 
then 

CPi — CP U {critical pair [V, Q] | Q G B} 
end if 

B< — B U {V} # no matter whether polj(V) ^ or not 
end if 
end while 

return {polynomial part of Q | Q G B} 
end 



Next is the main function of Reduction in ( Faugere . 2002| ) (function [2]). 
We can also write F5-reduction as a function which is similar to the TopReduction (func- 
tion EJ. 

Thus, the F5-reduction can be used to replace the TopReduction function perfectly in 
the Reduction function by simply updating the following codes: 

(Done', Todo')i — F5- Reduction^, B). 

Next, let us compare TopReduction and F5-reduction. Let Q G B, v = lm(J r )/lm(^), 
and suppose vQ is neither divisible nor rewritable by B or R. 

(1) For the case sign(J-") > sign(t> Q), there is no difference between the TopReduction and 
F5-reduction. 

(2) For the case sign(J r ) = sign(t>C?): 

• In TopReduction, labeled polynomial vQ is rewritable by B U {J 7 } . Thus the la- 
beled polynomial Q can never be selected from B, and nothing is done about Q in 
TopReduction. 



S 



Function 1 — TopReduction 



Input: a labeled polynomial J 7 , a set of labeled polynomials B, a set of labeled 
polynomials Todo which needs to be reduced further. 

Output: a 2-tuple (Done' , Todo'): both Done' and Todo' are sets of labeled poly- 
nomials, while Done' contains the result of reduction and Todo' includes 
the labeled polynomials to be reduced further. 

begin 

Qi — a labeled polynomial in B such that: 

(1) lpp((?) | lppOF) and denote v = lm(.F)/lm(£), 

(2) vQ is not divisible by B, and 

(3) vQ is not rewritable by B U {J 7 } U Todo. 
if such Q does not exist, 

then 

return ({J 7 }, 0) 
else 

if signature sign(J r ) >sign(i>C?), 
then 

return (0, {J 7 — t>£)} 
else 

return (0, {J 7 , vQ - J 7 )} 
end if 
end if 
end 



• In F5-reduction, no matter which labeled polynomial is used to F5-reduce J 7 , the 
reduction result of J 7 has the same signature with J 7 and will add to the set B. So when 
discussing the critical pair [J 7 , Q\ in later loops, this critical pair meets the Rewritten 
Criterion and is hence rejected. Thus no new labeled polynomial is generated. 

(3) For the case sign(J 7 ) < fsign(^): 

• In TopReduction, the labeled polynomial v Q — J 7 is calculated immediately and added 
to set Todo. Both labeled polynomials J 7 and vQ will be reduced in later loops of 
the Reduction function. Notice that the labeled polynomial J 7 may be selected from 
set Todo in later loops, but in that time, the labeled polynomial Q is not qualified to 
reduce J 7 , as v Q — T has been added to the set Todo and hence vQ is rewritable by 
B U {J 7 } U Todo. 

• While in F5-reduction, we cannot use the labeled polynomial Q to F5-reduce J 7 im- 
mediately. However, in the F5B algorithm, the critical pair [J 7 , Q\ must be added to 
the set CP, since labeled polynomial Q is in the set B. Moreover, the S-polynomial 
of the critical pair [J 7 , Q\ is exactly the labeled polynomial vQ — J 7 . The S-polynomial 
spoly(J 7 , Q) which can also be computed immediately in the next step if we use an 
appropriate strategy which selects the critical pair [J 7 , Q\ first. 
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Function 2 — Reduction 



Input: a set of labeled polynomials Todo which need to be reduced, a set of 

labeled polynomials B. 
Output: the reduction results of Todo reducing by the set B. 

begin 

Donei — 

while set Todo is not empty, 
do 

F< — the labeled polynomial with minimal signature in set Todo 
Todoi — Todo \ {J 7 }, 

(Done 1 , Todo')i — Top Reduction^, B, Todo) 

Donei — Done U Done' 

Todoi — Todo U Todo' 
end while 
return Done 
end 



Function 3 — F5-Reduction 

Input: a labeled polynomial J 7 , a set of labeled polynomisl B. 
Output: a 2-tuple (Done, Todo): both Done and Todo are sets of lebeled polyno- 
mials, while Done contains the result of reduction and Todo includes the 
labeled polynomials to be reduced further. 

begin 

Qi — a labeled polynomial in B such that: 

(1) lpp(^) | lpp(J r ) and denote v = lm(J r )/lm(^), 

(2) signature sign(J r ) >sign(vQ), 

(3) v Q is not divisible by B, and 

(4) vQ is not rewritable by B. 
if such Q does not exist, 

then 

return ({J 7 }, 0) 
else 

return (0, {T - vQ} 
end if 
end 



The third case sign(J r ) < vsiga(Q) is most complicated, so we illustrate the last case by 
a simple example. 

Example 4.2. Consider a system in Q[x, y, z] with the Graded Reverse Lex Order (x >- y >- 
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z): 

h = xz 2 + y 2 , 
f 2 = xy + xz, 
f 3 = yz + z. 



The corresponding labeled polynomials of fx, / 2 and /3 are Fx = (ei,/i), T% = (e 2 ,/ 2 ) 
and F 3 = (e 3 ,/ 3 ), respectively. The labeled polynomial set is B = {J 7 !, F2, F 3 }. We 
need to consider the critical pairs \Fx,F2\ = (v, Fx, z 2 , F2), [Fx,F 3 ] = (y,Fx,xz,F 3 ) and 
[F 2 , J3] = F 2 ,x,F 3 ) . 

Since the degree of lcm(lpp(J r 2), lpp(J"3)) = xyz is 3, which is lower than the other 
lcm's, we operate the critical pair [F2,F 3 ] first. Now we can see the major difference be- 
tween TopReduction and F5-reduction, since there exists a labeled polynomial Fx such that 
lpp(J-i) = xz 2 I xz 2 = lpp('P) but signature sign('P) <sign(7i). 



• In TopReduction, the labeled polynomial Q — Fx — P — [ex — ze 2 + a;e 3 , 22: + y 2 ) 
is calculated. Next, both labeled polynomials P and Q add to the set Todo for 
further reductions. In the next loop of the Reduction function, the labeled poly- 
nomial P = (ze2 — xe 3 ,xz 2 — xz) is irreducible, since the labeled polynomial T\ 
is rewritable by R this time. The labeled polynomial Q = {b\,xz + y 2 ) is irre- 
ducible as well. Then both V and Q add to set Done and create new critical pairs 

P], [P, T2I [V, 7 3 ], [J 7 !, Q], [Q, 7 2 ], [Q, 73], [Q, P]}. Combined with previous un- 
operated critical pairs {[7i, TV], [7i,7 3 ]}, the F5 algorithm continues to operate on 
these critical pairs. 

• While in F5-reduction, the labeled polynomial P is not F5-reducible by set B = 
{7i, J-2, J-3}. So the labeled polynomial P adds to set B immediately, and creates 
new critical pairs {[7i, P], [P, 7 2 ], [P, 7 3 ]}. Notice that there are still two un-operated 
critical pairs {[7i,7 2 ], [7i,7 3 ]}. When using a strategy of selecting critical pairs that 
selects the critical pair [Fi^P] first, then the S-polynomial spoly(J-i, "P) = T\ — P is 
calculated, and obtain Q = T\ — P = (ex — ze 2 +xe 3 , xz+y 2 ), which is not F5-reducible 
by set {7i, 7 2 , 7 3 , P}. Then the new critical pairs {[Fx, Q], [Q, F 2 ], [Q, F 3 ], [Q, P]} are 
created as well. 

In both cases, the remaining critical pairs are {[P, F2], [P, F3), [Fx, Q], [Q, F2), [Q,F 3 ], 
[Q,P], [Fx, F2], [Fx, F 3 ]}, since critical pair [Fx,P] will be rejected by Rewritten Criterion 
in TopReduction. 

After all, TopReduciton is equivalent to F5-reduction and hence the F5 algorithm is 
equivalent to the F5B algorithm with an appropriate strategy of selecting critical pairs. 
However, the F5B algorithm is simpler and easy to understand and analyze. Moreover, 
by using different strategies (of selecting critical pairs), the F5B algorithm becomes the 
Incremental F5 algorithm, Non-incremental F5 algorithm and Matrix F5 algorithm which 
are introduced in the next section. 

Since the F5B algorithm is equivalent to the F5 algorithm but m uch simpler to analyze , 



we propose a new proof for the correctness of the F5B algorithm in (ISun and Wangi . 120101 ) 



This new proof does not depend on the strategies of selecting critical pairs, so it also proves 
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the correctness for the other versions of the F5 algorithm. From this new proof, we find 
that the key of the F5B or F5 algorithm is the special reduction procedure which is a one- 
way reduction, i.e. only labeled polynomials with lower signatures can be used to reduce 
labeled polynomials with higher signatures. This fact will be much clearer in the Matrix F5 
algorithm introduced later. 



5. Some Versions of the F5 Algorithm 

In this section, we discuss three versions of the F5 algorithm. They are Incremental F5 
algorithm, Non-incremental F5 algorithm and Matrix F5 algorithm. 

5.1. Incremental F 5 algorithm 

Incremental F5 algorithm is the original F5 algorithm presented in (iFaugerd . 120021 ) . If 

the F5B algorithm uses an appropriate strategy of selecting critical pairs, then the F5B 
algorithm becomes the Incremental F5 algorithm. 

To describe this strategy appropriately, we need some new notations. For a labeled 
polynomial J 7 = (f , /) whose signature is sign(J-") = lpp(f ) = x a ei, we define the index 
of T to be index (J 7 ) := %. Given a critical pair \T , Q], we also define indexQj 7 , Q]) : = 
min{index(J-'), index(^)}. Then the strategy for incremental F5B algorithms can be de- 
scribed as follows: 

j<— max{index([.F,0]) | [F,Q\ E CP}, 

cp< — select a critical pair from CP with index j. 

Let us see an easy example. Let / 2 , fz, f$} C K[X] be the initial polynomials. Then 
at the beginning, set CP contains all the critical pairs: [J-'i, J-2], [^1,^3}, • • [J"3, J"4]- By 
the above strategy, j = 3 at the beginning. So the critical pairs with index 1 or 2 (such 
as [J-2, J-4]) cannot be selected unless all the critical pairs with index 3 have been operated. 
This means the critical pairs with index 1 or 2 can only be operated when the Grobner basis 
of ideal (/ 3 , / 4 ) is obtained. Similarly, critical pairs with index 1 (such as [J-'i, T3)) can be 
seleceted/operated only when the Grobner basis of ideal (f2, f3, fi) is computed. 

5.2. Non-incremental F5 algorithm 

Faugere also presents a Non-incremental F5 algorithm in a recent conference (INSCRPT 
2008). For example, it can use a strategy of selecting critical pairs like: 

G^min{deg(lcm(lpp(^),lpp(£))) 1 [?,Q] e CP}, 

cp< — select a critical pair from CP with degree d. 

In this case, for example, the critical pair [J-x, J-3] with index 1 may be operated earlier 
(not necessarily after the Grobner basis of ideal (f2, ■ ■ ■ , f m ) is obtained). 

However, this Non-incremental F5 algorithm is not really non-incremental. Since all the 
critical pairs should be operated sooner or later, the output of algorithm still contains the 
Grobner bases of the ideal (/j, • ■ ■ , f m ) where 1 < i < m. To transform the F5 algorithm to 
a r eal non-incremental algorithm, we need to chan ge the order of signatures. This is detailed 
in f Ars and Hashemi) . 120091 : ISun and Wang) . 12010 ). 
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5.3. Matrix F5 algorithm 

The most efficient version of the F5 algorithm is the Matrix F5 algorithm, which is 
an F5 algorithm that utilizes the matrix techniques introduced by the F4 algorithm when 
reducing S-polynomials. The Matrix F5 algorithm improves the original F5 algorithm in data 
structure, so the Matrix F5 algorithm differs from the F5B algorithm only by a strat egy of 



selecting critical pairs as well. Let us see the following example from (iFaugerd . 120021 ) . 



Example 5.1. Consider the homogeneous system in ¥ 2S [x,y, z] with the Graded Reverse 
Lex Order (x >- y >~ z) by the Matrix F5 algorithm^ 



f 3 = x 2 + 18xy + 19y 2 + 8xz + 5yz + 7z 2 , 
f 2 = 3x 2 + 7xy + 22xz + Uyz + 22z 2 + 8y 2 , 
f x = Q X 2 + 12xy + Ay 2 + \Axz + 9yz + 7z 2 . 



In order to compute the Grobner basis of (fi, f 2 , / 3 ), we set T\ = (ei, fi), T 2 = (e 2 , f 2 ) 
and J-3 = (e 3 ,/ 3 ). Next consider the critical pairs [J 7 !,^] = (1, J~i, 1, ^2), [^l,^] = 
(1, J 7 !, 1, J-3) and [J-2, ^3] = (1, I; Like the F4 algorithm, we use the part 1 x / b 1 x 
/2j 1 x fz to build the matrix of degree 2 in order to reduce the S-polynomials generated 
from these three critical pairs together: 

x 2 xy y 2 xz yz z 2 
f 3 /l 18 19 8 5 7 
A 2 = f 2 3 7 8 22 11 22 
U \ 6 12 4 14 9 7 

and after triangulation of the matrix A 2 : 




h 

B 2 = fa 
fi 

and two "new" polynomials appear: / 4 = xy + 4yz + 2xz + 3y 2 — z 2 (J^ = (— e 2 + 3e 3 ,/ 4 ) and 
sign(J r 4 ) = e 2 ) and f 5 = y 2 - llxz - 3yz - 5z 2 (J* 5 = (-4e a + 16e 2 - e 3 , f 5 ) and sign(J r 5 ) = 
ei), which are the reduction results of the S-polynomials spoly(J r i, T2), spoly(J r i, J-" 3 ) and 
spoly(7" 2 , Ta). 

For simplifying the statement, we will use the signatures of polynomials to replace the 
signatures of the corresponding labeled polynomials without confusions in the following of 
this example. 

Notice that the signature of the polynomial / 4 is e 2 , which corresponds to the label on 
the left of that row (underlined f 2 in the matrix B 2 ). 

Also we remark that the underlined 18 is not reduced by / 4 since the signature of / 3 
is e 3 which is smaller than e 2 (the signature of / 4 ). While the underlined is reduced, 



4 Here we set parameter b = directly. 
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since ei >- e 2 . This shows that the reduction procedure in the F5 algorithm is a one-way 
reduction. 

The next step is to consider the newly generated critical pairs: [J-i, J-4] = (y, T\, x, J-4), 

[^4,^2] = {x,T±,y,T 2 ), [T^T 3 \ = {x,T±,y,T 3 ), [^*5,^4] = (x,J"* 5 ,y, J"* 4 ), [7" 5 , = (x 2 , 7*5, 

y 2 , Jl), [J^JS] = [x 2 ,T^y 2 ,T 2 ) and [J"* 5 , J" 3 ] = (x 2 , J"" 5 , y 2 , J"" 3 ). We select these pairs by 
degree and build the matrix A 3 of degree 3 in order to operate the following critical pairs 

[J 7 !, J4] = (y, J4), [^4,^2] = (a;, F^y, ^2), [^4,^3] = (a:, Tk, y, ^3), 

[^5,^*4] = (X, 7*5, y, J"" 4 ) 

together. We only need to consider the parts y x / 3 , y x / 4 , x x / 4 ,x x / 5 , since the parts 
y x / 2 ,y x /1 are rewritable by J-4 and J-5 respectively. 

Like the F4 algorithm, the parts y x f 3 ,y x f±,x x f±,x x f$ are the rows to be reduced 
in the matrix, and we also need to select rows that are use to reduce these rows. Since 
power products y 3 ,x 2 z,xyz,y 2 z appear in the parts yf 3 , yf^, x/ 4 , xf 5 , we should add parts 
yf 5 , xf 3 , zfi, zf 5 to matrix A 3 in order to eliminate these power products. 

Now we have the matrix A 3 of degree 3. This matrix is ordered by the signatures of each 
row, which are listed in the round brackets: 







x 2 


y xy 2 


y 3 


x 2 z 


xyz 


2 2 
y z xz 


yz 2 


zfs 


(*e 3 ) 




( 








1 


18 


19 8 5 


7 \ 


yh 


(2/ e 3) 




1 


18 


19 





8 


5 7 





2/4 


(ze 2 ) 
















1 


3 2 4 


22 


A 3 = yf4 


(j/e 2 ) 







1 


3 





2 


4 22 





x/ 4 


(xe 2 ) 




1 


3 





2 


4 


22 





zh 


(zei) 



















1 12 20 


18 


yh 


(yei) 










1 





12 


20 18 





xh 


(xei) 






1 





12 


20 


18 


J 



and after triangulation (ordered by the leading power products of each row): 



x 2 y xy 


2 


3 2 

y x z 


xyz 


y 2 z 


2 

xz 


yz 2 


z 3 


yh (ye 3 ) 


/ 1 


18 


19 


8 


5 


7 


\ 




yh (ye 2 ) 





1 


3 


2 


4 


22 







xh (xe 2 ) 








1 





8 1 


18 


15 




B 3 = zf 3 (ze 3 ) 








1 


18 


19 8 


5 


7 




zfi (ze 2 ) 











1 


3 2 


4 


22 




zh (zei) 














1 12 


20 


18 




yh (yd) 














1 


11 


13 




xf 5 (^ei) 
















1 


18 ) 





and the polynomials h — y 3 + 8y z + xz 2 + 18yz 2 + 15z 3 (J"* 6 = ((15x + 18?/ + I2z)e 2 + 
(x + 7y + 17*)e 3 , / 6 ) and sign(.F 6 ) = xe 2 ), f 7 = xz 2 + Wyz 2 + 13^ 3 (T 7 = {(lSy + I8z)e 1 + 
(10x + 9y + 20^)e 2 + (16x + 13y + 13z)e 3 , f 7 ) and sign(JV) = ye x ) and / 8 = y^ 2 + I82 3 
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(Ts = ((lis + lly - 3z)ei + (21y + 9z)e 2 + (3x + 7y - 3z)e 3 , f 8 ) and sign(J r 8 ) = xei) 
are the reduction results of the S-polynomials of degree 3. Notice that although lpp(/e) = 
y 3 = ylpp(/s), labeled polynomial T§ is not F5-reducible by T*>. Thus j% is still a "new" 
polynomial. 

Now the Rewritten Criterion is much clearer. When building the matrix A3, we list the 
signatures of each row in round brackets. Labeled polynomials with the same signatures 
will play the same role in the matrix, so among the labeled polynomials with the same 
signatures, it suffices to deal with the latest results (that is why we care about the creating 
order of labeled polynomials) |f| 

Also the one-way reduction is evident in the matrix B 3 . Let us see the row x/4 (xe 2 ). The 
underlined 0, are reduced by rows zf 3 (ze 3 ) and z/4 (ze 2 ) respectively, while underlined 
8, 1, 18 are not eliminated by rows zf$ (zei), yf 5 (yei) and xf$ (xei). The reason lies in 
the one-way reduction. More specifically, the signatures of rows zf 3 (ze 3 ) and zf^ (ze 2 ) are 
ze% and ze 2 , both of which are smaller than the signature xe 2 of row x/4 (xe 2 ). Thus, the 
rows zf 3 (ze 3 ) and z/4 (ze 2 ) are able to reduce row x/ 4 (xe 2 ). However, we have signatures 
zei, yei, xe x y xe 2 , so the rows zf 5 (zei), yf$ (yex) and xf 5 (xei) are not qualified to reduce 
the row 2/4 (xe 2 ). Remark that, since only the rows x/4 (xe 2 ), yf 5 (ye{) and xfs (xe\) are 
worth saving, the others rows are not fully reduced in matrix S3. 

However, we must realize that, although the two new criteria of the F5 algorithm could 
reject almost all useless computations, the one-way reduction results in a poorer efficiency 
of eliminating matrix than the F4 algorithm. So it is really difficult to tell which of the F4 
and F5 algorithms is faster the other, especially in large examples. 

6. Conclusions 

In this paper, we rewrite F5 as F5B in Buchberger's style. We show that the F5B 
algorithm is equivalent to the original F5 algorithm. It is very easy to understand and 
implement F5B. The key of the F5B algorithm is the F5-reduction, which is a one-way 
reduction according to the signatures and their generating orders. F5B is also a useful 
tool to analyze F5 and F5-like algorithms. Although the F5 algorithm has good theoretical 
results for avoiding useless computations, the one-way reduction slows down the efficiency 
of algorithm more or less. It is desirable to have a more careful analysis of this issue, and 
we hope to work on it in the future. 
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5 If we add the row x/3 (xe^) to matrix A3, replacing the polynomials 2/4, y/4, 21/4, z/5, y/5, z/5 by 
polynomials 2/2, y/2, x/2, zfi, yfi, xf\ will lead to the same triangular form B3. 
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